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The functional method to derive the fractional Fokker-Planck equation 
for probability distribution from the Langevin equation -with Levy stable 
noise is proposed. For the Cauchy stable noise -we obtain the exact station- 
ary probability density function of Levy flights in different smooth potential 
profiles. We find confinement of the particle in the superdiffusion motion 
-with a bimodal stationary distribution for all the anharmonic symmet- 
ric monostable potentials investigated. The stationary probability density 
functions sho-w po-wer-law tails, which ensure finiteness of the variance. By 
revie-wing recent results on these statistical characteristics, the peculiarities 
of Levy flights in comparison -with ordinary Bro-wnian motion are discussed. 
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1. Introduction 

Anomalous diffusion in the form of Levy flights appears in many phys- 
ical, chemical, biological, and financial systems Levy flights are 
stochastic processes characterized by the occurrence of extremely long jumps. 
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The length of these jumps is distributed according to a Levy stable statistics 
with a power law tail and divergence of the second moment. This peculiar 
property strongly contradicts the ordinary Brownian motion, for which all 
moments of the particle coordinate are finite. The presence of anomalous 
diffusion can be explained as a deviation of real statistics of fluctuations 
from Gaussian law, that it has lead to the generalization of the central limit 
theorem by Levy and Gnedenko [5]-[7|. The divergence of the variance 
of Levy flights poses some problems as regards to the physical meaning of 
these processes. However, recently the relevance of Levy motions appeared 
in many physical, natural and social complex systems. The Levy type statis- 
tics, in fact, is observed in various scientific areas. Among many interesting 
examples we cite here the subrecoil laser cooling [Sj-pUj. the diffusion by 
flows in porous media [IT], the fluctuations in plasmas [12], the molecular 
collisions [13], the spatial gazing patterns of bacteria [13], the flights of an 
albatross [TB], the long paleoclimatic time series of the Greenland ice core 
measurements [16) . and the financial time series |17j-|19|. Experimental 
evidence of Levy processes was also observed in the motion of single ion 
in a one-dimensional optical lattice ^0] and in the particle evolution along 
polymer chains [2H [22] . 

The problem of the barrier crossing in a bistable potential, the particle 
escape from a metastable state, and the first passage time density have 
been analyzed, recently, for Levy flights The main focus in these 

papers is to understand how the barrier crossing behavior, according to the 
Kramers law [36], is modified by the presence of the Levy noise. 

Levy flights are a special class of Markovian processes, therefore the 
powerful methods of the Markovian analysis are in force in this case. We 
mean a possibility to investigate the stationary probability distributions of 
super diffusion, the first passage time and the residence time characteristics, 
the spectral characteristics of stationary motion, etc. Of course, this type 
of diffusion has a lot of peculiarities different from those observed in normal 
Brownian motion. The main difference from ordinary diffusion consists 
in replacing the white Gaussian noise source in the underlieing Langevin 
equation with a Levy stable noise. 

In this paper we use functional approach to derive the Fokker-Planck 
equation, with fractional space derivative, directly from Langevin equation 
with a Levy stable noise source. Starting from this equation we find the ex- 
act stationary probability distribution (SPD) of fast diffusion in symmetric 
smooth monostable potentials for the case of Cauchy stable noise. Specif- 
ically, we consider symmetric potential profiles U{x) = 7x^™'/(2m) (with 
odd m = 2n + 1 and even m = 2n), describing the dynamics of overdamped 
anharmonic oscillator driven by Levy noise. We find that for Levy flights in 
steep potential well, with steepness greater or equal to four, the variance of 



Duhkov' 11-05-07 printed on October 5, 2008 



3 



the particle coordinate is finite. This gives rise to a confined superdiffused 
motion, characterized by a bimodal stationary probability density, as previ- 
ously reported in Refs. [37]-[39]- However, in previous Ref. [^9] the authors 
analyzed the properties of stationary probability distribution for nonlinear 
Levy oscillators and its bimodality as a function of the Levy index a, by 
finding that it is more pronounced for a = 1 (Cauchy stable noise) and 
becomes the Boltzmann stationary distribution in the limit of a — > 2. Here 
we analyze the SPDs as a function of a dimensionless parameter (3, which is 
the ratio between the noise intensity D and the steepness of the potential 
profile 7. We find that the SPDs remain bimodal with increasing (3 param- 
eter, that is with decreasing the steepness 7 of the potential profile, or by 
increasing the noise intensity D. 



2. Functional method to derive the fractional Fokker-Planck 
equation from Langevin equation with Levy stable noise 

Ditlevsen and Yanovsky with co-authors for the first time obtained the 
fractional Fokker-Planck equation directly from Langevin equation, by re- 
placing the white Gaussian noise with Levy stable noise HO] (see also |4H 
142]). However, some attempts were undertaken before in Ref. \43>\ I44j . 

The theory of Levy processes is closely linked to that of infinitely divis- 
ible distributions [45j-[47j. Therefore, starting from this link, we have re- 
cently developed a more general approach, based on the theory of infinitely 
divisible distributions and functional analysis, to derive the generalized Kol- 
mogorov equation for arbitrary non-Gaussian white noise source [48j. Here 
we obtain the equation for probability distribution from Langevin equation 
with Levy stable noise, by a different approach with respect to that reported 
in Ref. [lo]. 

Let us consider the anomalous overdamped motion in the potential pro- 
file U (x) 

^^=-U'{x)+L{t). (1) 

Here x{t) is the displacement of particle and L (t) is the symmetric a-stable 
Levy noise with the characteristic function of increments 

{exTp {ik [VL (t + At) -TIL (t)]}) = /^expi^ik J L(T)dr|^ 

= exp{-D|A;|"At}, (2) 

where riL{t) is a generalized Wiener process [481 H9] which derivative is the 
Levy stable noise {drj^/dt = L[t)). Here a is the Levy exponent (0 < a < 2) 
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and D is the intensity of Levy noise. The case a = 1 corresponds to Levy 
noise L (t) with symmetric Cauchy distribution. First of all, we calculate 
the characteristic functional of the noise L (t). 

According to the definitions of the characteristic functional of the ran- 
dom process L{t) and the Stiltjes integral we have 



t 



&t[u] = ^exp |zy"n (r) L (r) (ir|y = ^exp l^y" u (r) dryi (r) 
= ^exp |i J2 ^ (^fc) ('^fc) ~ (Tfc-i)] I ^ 

= i™n ( n (^^) ^'^k) - VL (Tfc-l)]} ) , 

^^^^ \k=l I 

where "Qk is some internal point of the time interval (T,fc„i, r^) , S-r = max Ar^, 

k 

Arfc = Tfc — Tk~i (tq = 0, r„ = t). Taking into account that the increments 
of non-overlapping time intervals of the generalized Wiener process 
are statistically independent and using Eq. ^ we obtain 

n 

Qt M = lim W (exp {iu {'dk) [vl (^fc) - Vl (Tfc-i)]}) 

n 

= hm l[exp{-D\u{^kTATk} (3) 

= exp|-i:>^lim^^jM(??fc)|" Arfcj . 
By using the definition of the Riemann integral we finally get 

et[n]=exp|-L» I |n(r)rdr|. (4) 

To derive the fractional Fokker-Planck equation from Langevin equa- 
tion ([1]) we need the functional correlational formula for symmetric a-stable 
Levy noise L (t). We start from the generalization of Furutsu-Novikov for- 
mula |50il51j for arbitrary non-Gaussian random process £,{t), obtained pre- 
viously in [52] , 



iu (t) 



{Rt[C + z]), (5) 



._S_ 

iSz 
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where Rt [^] is a functional of noise ^ (t), defined on the observation interval 
(0,t), z {t) is a deterministic function, and $f [u] = ln6( [u]. Following Kly- 
atskin, we use the translation functional operator and taking into account 
that the function z{t) is deterministic we have 

(t) Rt[C + z]) = it) exp || e (r) ^^^1 ) [z] . (6) 
For the average entering in Eq. ([6]), after evident rearrangements, we find 



iu (t) dt 

Qt [u] d 
iu (t) dt 

[u] ' 
iu (t) 



In et[n] (7) 



exp < i J ^{t)u (r) dr 



Substituting Eq. ([7]) in Eq. and using again the functional translation 
formula we obtain Klyatskin result ([5]). By using the following integral 
representation for 

+ 00 

_ r (a + 1) sin (7ra/2) r 1 - cos(xm) 
m - - J -0+^ ax, [H) 



— oo 



we rewrite Eq. ^ as 

t +00 



e,[u] = e^p{-Qj dr J L^^(^^dx}, (9) 
where 

^ ^ DT (a + l)sin(7ra/2) ^ ^^^^ 
vr 

We obtain, therefore, the following expression for the variational operator 
in Eq. ([5]) for Levy stable noise L{t) 

*''"'-e/^=:^^-Q/^/e»'"»*. (11) 



iu{t) J iu{t)\x 



oo 
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Substituting this expression in Eq. ([5]) we get 

+00 X 

{L{t)Rt[L + z]) = Q J Jexp!^yj^yRt[L + z])dy. (12) 

—00 

By inserting the operator of functional differentiation into the average and 
by putting z = 0, we get finally 

+00 X 

{L (t) Rt [L]) =Q J -p^ J ^exp [vj^] Rt [L]) dy. (13) 
-00 

Now we are ready to derive the fractional Fokker-Planck equation for 
Levy flights using the functional approach. By differentiating, with respect 
to time t, the expression for probability density of random process x{t) 

W{x,t) = {6{x-x{t))) , (14) 

and taking into account Eq. dTJ, we obtain 

dW d d 

— = - [U'{x)W] --{L{t)5{x-x m . (15) 
To evaluate the average in Eq. (fTSi) we apply the formula (fT3l) 



{L{t)6ix-xm = Q I ^^l(exp!^yj^'^5{x-x{t))jdy. 

(16) 

Using functional differentiation rules, from Eq. ([1]) we get 

8{x-x {t)) = -—8{x-x (t)) = - 5(x-x (t)) . (17) 



5L{t) ' ^ " dx ' ^ " 6L{t) dx 

Thus, the variational operator 5/5L (t) with respect to the functional 5 {x — x (t)) 
is equivalent to the ordinary differential operator —d/dx. As a result, we 
have 

+00 z 

{L{t)5{x-x{t))) = Q j Jexpl^-y-^'jdyW{x,t). (18) 



-00 



After substitution of Eq. (jlSp in Eq. (jlSp and evaluation of the internal 
integral we arrive at 



d 



dz 

Wix,t)-^. (19) 

\z\ 
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By using the property of the translation operator 

exp|-z^|/(x) = /(x-z), (20) 

we arrive at the fohowing Kohnogorov equation for the probabihty density 
of nonhnear systems ([1]) driven by a symmetric a-stable Levy noise 

'-^-lW'im'].Qj ^^^--p''^''^ ^.. (21, 

— oo 

The Eq. (j2ip represents the well-known Fokker-Planck equation with frac- 
tional space derivative, which describes super diffusion in the form of Levy 
flights 

= 1 (22, 

3. Stationary probability distributions for Levy flights 

First of all, we can try to evaluate the stationary probability distribution 
Wst (x) from Eq. (j22p . if it exists. Of course, this evaluation is impossible for 
any potential profile, but the potential U{x) should satisfy some constraints. 
It is better to apply Fourier transform to the integro-differential equation 
(j2ip and to write the equation for the characteristic function 

+ 00 

^{k,t) = (e'^<'^ ) = J e'^'^'W {x, t) dx. (23) 

— oo 

After simple manipulations we find 

+ 00 

-g^ = -ik I e''''^U'{x)W{x,t)dx-D\k\'''&. (24) 

— oo 

For smooth potential profiles U (x), expanding in power series near the point 
X = 0, we can rewrite this equation in the operator form 

-ikU' (-i-^)^-D\k\"^. (25) 



dt V dk_ 

In particular, for stationary characteristic function, from Eq. (j25p we get 

U' (-i^) ^st - iD I fc^-i sgn A: • = , (26) 
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where sgnA; is the sign function. Unfortunately, one cannot solve Eq. (I26p 
for arbitrary potential U (x) and arbitrary Levy exponent a. 

Let us consider, as in [39], the symmetric smooth monostable potential 
U (x) = 7x^™/ (2m) (m = 1, 2, . . .). The Eq. ()26p . therefore transforms into 
the following differential equation of (2m — l)-order 

+ (-1)"+! /32™-i I A;r-i sgn • = , (27) 



where /3 = ^"^'-i/lTpy. As it was proved by analysis of Eq. ([27|) in [37] . 
the stationary probability distribution Wst (x) has non-unimodal shape and 
power tails 

Wstix)r^ .aJ+^-i . Ixl^oo. (28) 
|x| 

In Ref . [37] , the estimation of bifurcation time for transition from unimodal 
initial distribution to bimodal stationary one and the existence of a transient 
trimodal state for m > 2 were found. 

Exact solution of Eq. (I27p can be only obtained for the case of Cauchy 
noise: a = 1. Due to the symmetry of the characteristic function 'dgt (—k) = 
'&st (k) we can reduce Eq. ([27]) to the linear differential equation with con- 
stant parameters 



St (_i)"^^2m-i^^^^Q (A:>0). (29) 



From the corresponding characteristic equation 

we select the roots with negative real part, which are meaningful from phys- 
ical point of view. The general solution of Eq. (j29p . therefore, reads 

^st{k)= E vlzexp|-/?|fc|cos ^^^^_^ 

( . 7r(m-2/-l) \ 
cos(/3|A;|sm 2m - I '^V ' ^ ' 

where the quadratic brackets in the upper limit of the sum denote the integer 
part of the expression. The unknown constants Ai and (pi can be calculated 
from the obvious conditions 

T?,aO) = l, 4t'"'^(+0) = (j = l,2,...,m-l). (32) 
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Substituting Eq. ()3ip in Eq. ([32]) we have 



[(m-l)/2] 

^ Al cos = 1 

[(m-l)/2] 

X! COS 



7r(2j - 1) (m + 2/) 
2m - 1 



(33) 

(j = l,2,...,m-l). 



Making the reverse Fourier transform of Eq. (j3ip we obtain the stationary 
probability distribution (SPD) of the particle coordinate 



7r(m-2«-l) , 
2m-l + 



+ /J^ COS 



7r(m-2«-l) 

2m-l - '^l 



2x2/32 COS ^g±ii+/?4 



[(m-l)/2] ^2 

(34) 

The parabolic potential profile U {x) = 7^2/2 corresponds to a linear 
system ([T]). In this situation, from Eqs. ()33p and ()34p we easily obtain the 
following obvious result 



W,t ix) 



P 



7r(x2 + /?2) 



(35) 



i.e. due to the stability of the Cauchy distribution (|35l) . the probabilistic 
characteristics of driving noise increments (see Eq. ([2])) and Markovian pro- 
cess x (t) are similar. 

For quartic potential {m = 2), from the set of Eq. (j33p . we find = 
2/\/3, = 7'"/6. Substituting these parameters in Eq. (13^ we obtain 



W,t (x) 



vr (x4 - x2/32 + /?4) ' 



(36) 



which coincides, for /? = 1, with the result obtained in Ref. ^39j- The plots 
of stationary probability distributions ()36p for Levy flights in symmetric 
quartic potential for different values of parameter f3 are shown in Fig. 1. The 
superdiffusion in the form of Levy flight gives rise to a bimodal stationary 
probability distribution when the particle moves in a monostable potential, 
differently from the ordinary diffusion of the Brownian motion characterized 
by unimodal SPD. 

The SPD of superdiffusion has two maxima at the points x = ib/3/-v/2, 
with the value (VF^j)^^^ = 4/ (37r/?). Since the value of the minimum is 
Wst (0) = 1/ (vr/3), the ratio between maximum and minimum value is con- 
stant and equal to 4/3. The width of probability density increases with 
increasing parameter P = ^/ITpy, i. e. with decreasing the steepness 7 of the 
quartic potential profile, or with increasing the noise intensity D. 
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Fig. 1. Stationary probability distributions for Levy flights in symmetric quartic 
potential U{x) — 'yx^/A for different values of dimensionless parameter (3: 1 - 
P = 0.5, 2 - /3 = 1, 3 - /? = 1.5. 



Carrying out analogous procedure we obtain the stationary probability 
distributions for the cases m = 3, 4, 5 



Wst {x) 
Wst {x) 
Wst ix) 



(37) 



TT (x2 + /32) - 2P^X^ COS 7r/5 + p^) ' 

p[ 

TT (X^ - 2/32x2 cos 7r/7 + + 2/32x2 COS 27r/7 + (3^) ' 

^! 

TT (x2 + /32) (x4 - 2/32x2 COS 7r/9 + /34) (x4 + 2/^2x2 COS 47r/9 + ■ 



The plots of distributions (j37p . for different values of parameter /?, are re- 
spectively shown in Figs. 2-4. It must be emphasized that according to 
Figs. 2-4, these distributions remain bimodal and have the same tendency 
with increasing /?, but the ratio between maximum and minimum increases 
with increasing m. From Eqs. (j36p and (j37p we see that the second mo- 
ment of the particle coordinate is finite for m > 2. This means that there 
is a confinement of the particle motion due to the steep potential profile, 
even if the particle moves according to a superdiffusion in the form of Levy 
flights [37J. The presence of two maxima is a peculiarity of the superdiffu- 
sion motion. Because of the fast diffusion due to Levy flights, the particle 
reaches very quickly regions near the potential walls on the left or on the 
right with respect to the origin x = 0. Then the particle diffuses around 
this position, until a new flight moves it in the opposite direction to reach 
the other potential wall. As a result, the particle spends a large time in 
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Fig. 2. Stationary probability distributions for Levy flights in symmetric potential 
U{x) = 7x^/6 for different values of dimensionless parameter (3: 1 - (3 = 0.5, 2 - 
/3 = 1, 3 - /3 = 1.5. 



W,lx) 




-2-10 12 



Fig. 3. Stationary probability distributions for Levy flights in symmetric potential 
U{x) = for different values of dimensionless parameter (3: \ - (3 — 0.5, 2 - 

/3 = 1, 3 - /3 = 1.5. 

some symmetric areas with respect to the point x = 0, differently from the 
Brownian diffusion in monostable potential profiles. These symmetric areas 
lie near the maxima of the bimodal SPD. For fixed D and m, these maxima 
are closer or far away the point x = depending on the greater or smaller 
steepness 7 of the potential profile. This corresponds to a greater or smaller 
confinement of the particle motion. Of course, such a confinement is more 
pronounced for greater m, that is for steeper potential profiles. 

On the basis of Eqs. (|35p - (|37p and the known behavior of density tails 
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Fig. 4. Stationary probability distributions for Levy flights in symmetric potential 
U{x) = 7a;"'^*'/10 for different values of dimensionless parameter (3: 1 - 13 — 0.5, 2 
- /3 = 1, 3 - /3 = 1.5. 



), we can write the general expressions for stationary probability distri- 
bution in the case of potential U (x) = 73;^™/ (2m) with odd m = 2n + 1 



Qin+l n-1 

(x) = _ , n 



1 



vr (x2 + /?2) 11 - 2/32^2 cos [vr (4/ + 1) / (4n + 1)] + (3^ ' 



and even m = 2n 

o4n-l n-1 



Wst (x) 



n 



^ fj-^ - 2/32x2 COS [vr (4/ + 1) / (4n - 1)] + /^^ ' 
which are, together with Eqs. (j36p - (j37p . the main result of this paper. 



(38) 
(39) 



4. Conclusions 

We used functional analysis approach to derive the fractional Fokker- 
Planck equation directly from Langevin equation with symmetric a-stable 
Levy noise. This approach allows to describe anomalous diffusion in the 
form of Levy flights. We obtained the general formula for stationary prob- 
ability distribution of superdiffusion in symmetric smooth monostable po- 
tential for Cauchy driving noise. All distributions have bimodal shape and 
become more narrow with increasing steepness of the potential or with de- 
creasing noise intensity. We found that the variance of the particle coordi- 
nate is finite for quartic potential profile and for steeper potential profiles, 
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that is a confinement of the particle in a super diffusion motion in the form 
of Levy flights. As a result, we can evaluate the power spectral density of 
a stationary motion. Calculations of residence times for the case of Levy 
flights in bistable potential with steep potential wells, and anomalous dif- 
fusion in periodic ratchet-like potentials are the subjects of forthcoming 
investigations. 

This work has been supported by MIUR, CNISM, and by Russian Foun- 
dation for Basic Research (project 05-02-16405). 
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